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Abstract. Abundances of 22 chemical elements, including the key elements 
and isotopes such as ^^C, ^^C, N and O, are investigated in the spectrum of 
33 Psc, a single-lined RS CVn-type binary of low magnetic activity. The high 
resolution spectra were observed on the Nordic Optical Telescope and analyzed 
with the MARCS model atmospheres. The following main parameters have 
been determined: T^s = 4750 K, log g = 2.8, [Fe/H] = -0.09, [C/Fe] = -0.04, 
[N/Fe] = 0.23, [0/Fe] = 0.05, C/N = 2.14, i^c/^C = 30, which show the 
first-dredge-up mixing signatures and no extra-mixing. 

Key words: stars: RS CVn binaries, abundances - stars: individual (33 Psc 
= HD 28)) 

1. INTRODUCTION 

This is the third paper in a series dedicated to a detailed study of photospheric 
abundances in RS CVn stars (Tautvaisiene et al. 2010; Barisevicius et al. 2010, 
hereafter Papers I and II) with the main aim to get the carbon isotope ^^C/^'^C 
and C/N ratios in these chromospherically active stars. We plan to investigate 
correlations between the abundance alterations of chemical elements in the atmo- 
spheres of these stars and their physical macro parameters, such as the speed of 
rotation and the magnetic field. 

Here we present results of the analysis of a bright KOIIIb single-line binary 
33 Psc (HD 28) of = 4.78 mag (van Leeuwen 2007). As a binary with an orbital 
period of P = 72.93 days it was recognized by Harper (1926). The presence of 
the emission in its H and K lines of Ca ii was detected spectroscopically by Young 
& Konigest (1977). From the analysis of intensity of the Can emission in the K 
line, Glebocki & Stawikovski (1979) concluded that the primary is of 1.6 Mq and 
the secondary is of 0.7 Mq. This evaluation was confirmed by Pourbaix & Boffin 
(2003) using the Hipparcos intermediate astrometric data: the mass of the primary 
was found to be 1.7 ± O.4M0 and of the secondary 0.76 ± 0.11 AIq. This system 
has not been resolved by speckle interferometry (Hartkopf et al. 2001). 

From narrow-band photometry, Hansen & Kjaergaard (1971) and Gottlieb & 
Bell (1972) have determined that 33 Psc is slightly metal-deficient ([Fe/H] « -0.2). 
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The high-resolution spectroscopic metallicity determinations of 33 Psc are sUghtly 

controversial: McWiniam (1990) has obtained [Fe/H] w-0.31, while the [Fe/H] 
values obtained by Randich et al. (1994) and Morel et al. (2004) are solar. 

In papers I and II wc investigated A And and 29 Dra which are among most 
active RS CVn stars. 33 Psc by now has the lowest known magnetic activity among 
RS CVn stars (Strassmeier et al. 1988). Results of Einstein X-ray observations of 
33 Psc were presented by Walter (1985). Its surface flux is < 5 • 10'^ crgcm"^ s~^, 
which is about a factor of three lover than the flux from the quiet solar corona, and 
is two orders of magnitude below the surface fluxes observed in the least-active 
long-period RS CVn systems. Drake et al. (1989) for 33 Psc have found quite low 
upper limit of the 6 cm radio flux, implying low activity as well. The ultraviolet 
chromosferic emission lines are also weak (Basri et al. 1985). 

For A And and 29 Dra the ratios of ^^C/^^C were found to be lower than it is 
predicted by the first dredge-up theory of the first-ascent giants. Therefore it was 
interesting to check what are the effects of mixing processes in 33 Psc, which lies 
on the giant branch almost at the same luminosity but is much less active than 
A And and 29 Dra. 

2. OBSERVATIONS AND THE METHOD OF ANALYSIS 

The spectrum of 33 Psc was observed twice - in August of 1999 and Septem- 
ber of 2006 on the 2.56 m Nordic Optical Telescope using the SOFIN cchelle 
spectrograph with the optical camera, which provided a spectral resolving power 
of i? 80 000, for 26 slightly shifted spectral orders, each of ~ 4 nm, in the spec- 
tral region from 500 to 830 nm. Details of spectral reductions and the method of 
analysis are presented in Paper I. 
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Fig. 1. A comparison of the equivalent widths measured in the 33 Psc spectra 
observed in 1999 and 2006. 
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Table 1. Equivalent widths of lines, EW, in the spectrum of 33 Psc. 



Element 


A (A) 


EW{m.k) 


Element 


A (A) 


EW(rak) 


Element 


A (A) 


EW{m.k) 


Sil 


5517.55 


23.4 




6274.66 


74.0 




6392.54* 


62.3 




5645.60 


52.0 




6285.16 


83.1 




6574.21* 


93.1 




5793.08 


56.0 




6292.82 


87.1 




6581.21* 


75.9 




5948.54 


91.6 


Cn 


5712.78 


49.9 




6646.97 


43.0 




6131.85 


29.5 




5783.87 


76.6 




6786.86* 


48.8 




7003.57* 


51.7 




5784.97 


63.2 




6793.27* 


31.3 


Cai 


5260.38 


66.9 




5787.92 


81.5 




6839.83 


78.7 




5867.57 


53.1 




6661.08 


27.5 




6842.69 


60.4 




6455.60 


94.5 




6979.80 


73.6 




6843.65 


79.7 




6798.47* 


23.7 




6980.91 


27.0 




6851.64 


33.3 


Sen 


5526.81 


106.2 


Fel 


5395.22 


42.3 




6857.25 


45.9 




5667.14 


67.5 




5406.78 


58.9 




6858.15 


71.6 




6279.75 


61.2 




5522.45 


65.8 




6862.49 


51.0 




6300.69 


22.6 




5577.03 


21.1 




7461.53* 


71.7 


Til 


5648.58 


44.7 




5579.35 


23.8 


Fell 


5132.68 


33.5 




5662.16 


66.6 




5607.67* 


36.3 




5264.81 


47.0 




5716.45 


33.7 




5608.98* 


29.4 




6113.33 


16.0 




5739.48* 


39.6 




5651.48 


37.3 




6369.46 


23.6 




5880.27 


58.0 




5652.33 


48.4 




6456.39 


55.2 




5899.30* 


95.6 




5653.86 


57.6 




7711.72* 


43.7 




5903.31* 


44.7 




5679.03 


73.8 


Co I 


5647.23 


60.0 




5941.76 


81.3 




5720.90 


33.1 




6117.00 


40.8 




5953.17 


87.7 




5732.30 


23.9 




6188.98 


58.0 




5965.83 


83.1 




5738.24* 


32.2 




6455.00 


44.4 




6064.63* 


58.5 




5741.86* 


52.6 




6678.82 


42.4 




6098.66 


25.3 




5784.67 


57.8 


Nil 


5587.87* 


98.1 




6121.00 


27.6 




5793.92 


56.9 




5589.37* 


45.0 




6126.22 


83.8 




5806.73 


72.6 




5593.75* 


63.1 




6220.49* 


40.7 




5807.79 


30.7 




5643.09 


30.3 




6303.77 


52.4 




5811.92 


27.1 




5748.35* 


77.7 




6599.11* 


66.7 




5814.82 


45.7 




5805.22 


58.8 




6861.45 


41.9 




6027.06 


86.7 




6053.68* 


36.3 


Vi 


5604.96* 


39.2 




6034.04 


21.1 




6111.08 


54.9 




5646.11 


42.7 




6035.35 


17.4 




6128.98 


73.8 




5657.45 


52.8 




6054.07* 


23.5 




6130.14 


37.3 




5668.37 


52.8 




6056.01* 


86.8 




6204.60 


43.2 




5743.43* 


60.2 




6098.25 


32.1 




6378.25 


51.4 




6039.74* 


69.2 




6105.13 


24.0 




6598.60* 


41.8 




6058.18* 


35.2 




6120.24 


47.4 




6635.13 


44.9 




6111.65 


75.2 




6187.99 


75.6 




6767.78 


132.3 




6119.53 


82.9 




6200.32 


114.1 




6772.32* 


74.9 




6135.37 


69.8 




6226.74* 


54.2 




6842.03 


51.6 




6224.50* 


68.6 




6229.23* 


75.7 




7001.55* 


49.1 




6233.19* 


62.5 




6270.23 


92.5 




7062.97 


57.9 




6266.30 


47.9 




6380.75 


76.8 




7715.59* 


78.7 



★ The asterisk indicates the averaged equivalent widths from the spectra of 1999 and 2006, 
while other equivalent widths are only from the spectrum of 1999. 



In the spectra we selected 135 atomic lines for the measurement of equivalent 
widths and 17 lines for the comparison with synthetic spectra. The measured 
equivalent widths of lines are presented in Table 1, in overlapping regions they 
show a very good agreement (see Figure 1) and were averaged. 
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2.1. Atmospheric parameters 

Initially, the effective temperature, T^s, for 33 Psc was estimated and averagcid 
from the intrinsic color indices {B — V)o and (6 — y)o, using the calibrations by 
Alonso et al. (1999). The values of color indices, B-V = 1.06 and 6 - y = 0.627, 
were taken from Reglero et al. (1987) and Hauck & Mermilliod (1998), respectively. 
A small dereddening correction of Eb-v = 0.01, estimated using the Hakkila et 
al. (1997) software, was taken into account. 

The agreement between the temperatures deduced from the two color indices 
was quite good, with a difference of 80 K only. No obvious trend of the Fel 
abundances with the excitation potential was found (Figure 2). 

The surface gravity log g was found by adjusting the model gravity to yield the 
same iron abundance from the Fei and Fe ii lines. The microturbulent velocity 
Vt corresponding to a minimal line-to-line Fe i abundance scattering was chosen as 
the correct value. Consequently, [Fe/H] values do not depend on the equivalent 
widths of lines (Figure 3). 




Fig. 2. The [Fcl/H] abundance values versus the lower excitation potential xi- 
The mean abundance ([Fel/H] = -0.09 dcx) is shown as a dotted line. 




60 
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Fig. 3. The [Fei/H] abundance values versus the equivalent widths. The mean 
abundance ([Fei/H] = -0.09 dex) is shown as a dotted line. 
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2.2. Mass determination 

The mass of 33 Psc was evaluated from its effective temperature, luminosity 
and the isochrones from Girardi et al. (2000). The luminosity log(L/L0) = 1.39 
was calculated from the Hipparcos parallax tt = 25.32 mas and V = 4.61 mag 
(van Leeuwen 2007), the bolometric correction calculated according to Alonso et 
al. (1999), and the above mentioned Eb~v — 0.01. The mass of 33 Psc ^ 1.6 M© 
was found. 

Previous mass determinations for 33 Psc are within 1.1 (Gottlieb & Bell 
1972) and 3.0 M© (Barrado y Navascues et al. 1998). The masses 1.6 M© obtained 
by Glebocki & Stawikowski (1979) and 1.7 M© by Pourbaix & Boffin (2003) and 
Morel et al. (2004) are close to our result. 



2.3. Spectrum syntheses 

The method of synthetic spectra was used to determine carbon abundance 
from the C2 line at 5135.5 A. The interval 7980-8130 A, containing strong ^^c^N 
and ^^C-'^^N features, was used for the nitrogen abundance and the ^^C/^'^C ratio 
determinations. The ^^C/^'^C ratio was determined from the (2,0) ^''C^^N feature 
at 8004.7 A. All log gf values were calibrated to fit to the solar spectrum by Kurucz 
(2005) with solar abundances from Grevesse & Sauval (2000). 

The oxygen abundance was determined from the forbidden [O i] line at 6300.31 A 
with the oscillator strength values for ^®Ni and ^''Ni from Johansson et al. (2003) 
and the values log gf = —9.917 obtained by fitting to the solar spectrum (Kurucz 
2005) and log A© = 8.83 taken from Grevesse & Sauval (2000). 

In Figures 4 and 5 we show several examples of synthetic spectra in the vicinity 
of C2, [Oi] and i^ci^N lines. 
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Fig. 4. Synthetic spectrum fits to the C2 line at 5135.5 A and the forbidden [Ol] 
line at 6300.3 A. The observed spectra are shown as solid lines. The dashed lines are 
synthetic spectra with [C/Fe] = -0.03, 0.13 and 0.23 dex and [0/Fe] = 0.03, 0.13 and 
0.23 dex. 
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Fig. 6. Synthetic spectrum fit to the Ndll hne at 5276.806 A, Yll at 5402.78 A 
and La II at 6390.48 A. The observed spectrum is shown as a soUd hne. The dashed 
lines are synthetic spectra with [Nd/Fe] = 0.11, 0.21 and 0.31 dex, [Y/Fe] = -0.07, 
0.03 and 0.13 dex and [La/Fe] = -0.07, 0.03 and 0.13 dex for Ndii, Yii and Laii lines, 
respectively. 

The abundance of Nai was estimated using the line at 5148.84 A which is 
slightly blended by the Nil 5148.66 Aline. These two lines are distinct in the 
Sun, so we were able to calibrate their log gf values using the solar spectrum. 
So, the sodium abundance value in our study is affected by uncertainty of the 
nickel abundance determination by the Equivalent Widths method. Fortunately, 
the line-to- line scatter of [Ni/H] determination from 20 lines of Nil was as small 
as 0.06 dex. 

For the evaluation of Zri abundance the lines at 5385.13 A, 6127.48 A and 
6134.57 A were used. Evaluation of the Yii abundance was based on the 5402.78 A 
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line, La II on the 6390.48 A line, Cell on the 5274.22 A and 6043.38 A lines and 
Ndii on the 5276.86 A line. The synthetic spectrum fits to the Ndii line at 
5276.806 A, Yii at 5402.78 A and Laii at 6390.48 A arc displayed in Figure 6. 

The abundance of r-process element praseodymium was based on the Pr ii line 
at 5259.72 A and of europium on the Euii line at 6645.10 A (Figure 7). The 
hyperfine structure of Eu II was taken into account when calculating the synthetic 
spectrum. The wavelength, excitation energy and total log gf = 0.12 were taken 
from Lawler et al. (2001), the isotopic meteoritic fractions of ^^^Eu, 47.77%, and 
i^^Eu, 52.23%, and isotopic shifts were taken from Biehl (1976). 

Due to slow rotation spectral lines of 33 Psc are broadened very little. We 
used tisini = 1.9 kms^^ from Batten et al. (1989), which fits well the profiles of 
spectral lines of 33 Psc. 




5259.5 5260.0 6645.0 6645.5 



1(A) 

Fig. 7. Synthetic spectrum fit to the Euii line at 6645.10 A. The observed spec- 
trum is shown as a solid line. The dashed lines are synthetic spectra with [Eu/Fe] = 0.02, 
0.12 and 0.22 dex. 

2.4- Estimation of uncertainties 

The sources of uncertainty were described in detail in Paper I. The sensitivity 
of the abundance estimates to changes in the atmospheric parameters for the 
assumed errors (± 100 K for Teg, ±0.3 dex for log g and ±0.3 km s~^ for Vt) is 
illustrated in Tabic 2. It is seen that possible parameter errors do not affect the 
abundances seriously; the element-to-iron ratios, which we use in our discussion, 
are even less sensitive. 

The scatter of the deduced line abundances cr, presented in Table 3, gives 
an estimate of the uncertainty due to the random errors, e.g., in the continuum 
placement and the line parameters (the mean value of a is 0.06 dex). Thus, the 
uncertainties in the derived abundances originating from the random errors are 
close to this value. 

Since the abundances of C, N and O arc bound together by the molecular 
equilibrium, we have also investigated how the error in one of them typically 
affects the abundance determination of others. 

A[0/H] = 0.10 causes A[C/H] = 0.04 and A[N/H] = 0.05; A[C/H] = 0.10 
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causes A[N/H] = -0.12 and A[0/H] = 0.02. A[N/H] = 0.10 has no effect on both 
carbon and oxygen abundances. 

Table 2. Sensitivity of tfic abundances to changes of the atmospheric parameters. 
The table entries show the effects on the logarithmic abundance relative to hydrogen, 
A[E1/H]. 



Ion 




Alog g 




Ion 


AT.ff 


Alog g 


Aut 




+100 K 


+0.3 


+0.3km s^l 




+100 K 


+0.3 


+0.3km s-l 


C(C2) 


0.02 


0.11 


0.00 


Fell 


-0.09 


0.17 


-0.06 


N(CN) 


0.05 


0.13 


0.00 


Co I 


0.05 


0.06 


-0.08 


0([0i]) 


0.00 


0.13 


0.00 


Nil 


0.01 


0.06 


-0.10 


Nai 


0.08 


-0.01 


0.00 


Yii 


0.00 


0.11 


-0.03 


Si I 


-0.04 


0.06 


-0.04 


Zrl 


0.20 


-0.04 


-0.01 


Cai 


-0.09 


-0.01 


-0.08 


Lao 


0.01 


0.13 


-0.01 


Sen 


-0.02 


0.12 


-0.09 


Ceo 


0.00 


0.08 


0.00 


Til 


0.12 


0.01 


-0.08 


Prn 


0.01 


0.13 


-0.01 


Vi 


0.15 


0.01 


-0.10 


Ndn 


0.01 


0.13 


-0.01 


Cri 


0.08 


-0.01 


-0.08 


Eun 


-0.01 


0.13 


0.00 


Fel 


0.04 


0.04 


-0.08 










C/N 


-0.23 


-0.10 


0.10 


12C/13C 


-2 


1 


-2 



Table 3. Element abundances relative to hydrogen, [El/H]. ct is a standard deviation 
in the mean value due to the line-to-line scatter within the species. A'^ is the number 
of lines used for the abundance determination. 



Ion 


AT 


[El/H] 




Ion 


N 


[El/H] 


cr 


C(C2) 


1 


-0.13 




Fen 


6 


-0.08 


0.06 


N(CN)* 


4 


0.14 


0.02 


Co I 


5 


0.05 


0.07 


0([0i]) 


1 


0.04 




Nil* 


20 


0.00 


0.06 


Nal 


1 


-0.01 




Yii 


1 


-0.06 




Si I* 


6 


0.04 


0.06 


Zri 


3 


-0.21 


0.08 


Cai* 


4 


-0.02 


0.07 


La II* 


1 


-0.06 




Sen 


4 


0.13 


0.08 


Cell 


2 


-0.12 


0.07 


Til* 


18 


-0.01 


0.06 


Pril 


1 


-0.02 




Vi* 


16 


0.07 


0.07 


Ndii 


1 


0.12 




Cri 


7 


0.01 


0.05 


Euii 


1 


0.03 




Pel* 


49 


-0.09 


0.04 











* The asterisk indieates the elements, wliosc abundances were calculated from the 1999 and 
2006 spectra, while the others were calculated only from the 1999 spectrum. 



3. RESULTS AND DISCUSSION 

The following atmospheric parameters are obtained for 33 Psc: Tefj = 4750 K, 

\ogg = 2.8, vt = 1.1 kms-i, [Fe/H] = -0.09, [C/Fc] = -0.04, [N/Fe] = 0.23, 
[0/Fe] = 0.13, as well as the ratios C/N = 2.14 and i^C/i^C = 30. The element 
abundances [A/H] and the a values (the line-to-linc scatter) are listed in Table 3. 

In Figure 8 the obtained values of [El/Fe] for 33 Psc are compared with the 
results of other authors. Although McWilliam (1990) has used a high-resolution 
spectrum, his values of [El/Fc] arc rather scattered, probably due to small number 
of lines measured. The recent results of Morel et al. (2004) are in good agreement 
with ours. The element to iron ratios in 33 Psc are close to solar. 
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Fig. 8. Element abundances for 33 Psc, as determined in this work (filled diar- 
monds), by Morel ct al. (2004. circles) and by McWilliam (1990. triangles). 




The evolutionary sequences in the luminosity vs. effective temperature diagram 
by Girardi et al. (2000) show that 33 Psc with its luminosity log(L/_L0) = 1.39 is a 
first ascent giant lying slightly below the red giant sequence bump at log (L/Lq) = 
1.6 (Charbonnel & Lagarde 2010). According to the mentioned model of mixing, 
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carbon and nitrogen abundances in 33 Psc should be altered only by the first 

drcdgc-up. Our results confirm this prediction. In Figure 9 we compare C/N 
and ^^C/^^C ratios of 33 Psc with the theoretical models including extra-mixing. 
The model, called 'cool bottom processing' (CBP), was proposed by Boothroyd & 
Sackmann (1999) and the model, called 'thermohaline mixing' (TH), was proposed 
by Charbonnel & Lagarde (2010). The position of 33 Psc in Figure 9 indicates 
that its ratios of carbon isotopes and C/N have not been altered by extra mixing. 

Two more active RS CVn stars A And and 29 Dra, investigated by now in 
our program, gave a hint that extra-mixing processes may start acting in these 
low-mass chromosphcrically active fast rotating stars slightly earlier than at the 
bump of the red giant sequence in non-active stars. The star 33 Psc with negligible 
activity seems to be normal in this respect. 

The abundance of lithium is also very sensitive to mixing. During the first 
dredge-up, for a star of the mass of 33 Psc, the Li abundance drops to approx- 
imately logA(Li) = 1.44 (Charbonnel & Lagarde 2010). However, the avail- 
able determinations of Li abundance in 33 Psc show lower abundances: Brown 
et al. (1989) found logA(Li) = 0.8, Pallavicini et al. (1992) and Randich et 
al. (1994) found log^(Li) < 0.1, and Barrado y Navascues et al. (1998) found 
logyl(Li) = 0.41. These values are lower than predicted by the first dredge-up 
model, while the carbon and nitrogen abundances are in agreement with it. Thus, 
the lowered lithium abundance should be related to other mechanisms. 
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